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Abstract 

We provide a polynomial time algorithm for computing the universal Grobner basis of 
any polynomial ideal having a finite set of common zeros in fixed number of variables. One 
ingredient of our algorithm is an effective construction of the state polyhedron of any member 
of the Hilbert scheme Hilb^ of n-long d-variate ideals, enabled by introducing the Hilbert 
zonotope TC^ and showing that it simultaneously refines all state polyhcdra of ideals on Hilbfj . 

1 Introduction 

The universal Grobner basis of an ideal I in the algebra ¥[x] := ¥[x±, . . . ,Xd] of <i-variate poly- 
nomials over a field is the minimal set U{I) which is simultaneously a Grobner basis for / under 
every monomial order. A finite universal Grobner basis always exists and in a sense is the ul- 
timate generating set of / for algorithmic purposes. In particular, for ideals having a finite set 
of common zeros (variety) over the algebraic closure of F, a universal Grobner basis reduces the 
problem of computing the zero set to the problem of finding roots of d univariate polynomials. 
For instance, consider the system P = {x\ — xi-,x\ — 7x2 + 6xi,xiX2 — 3x2 + 2xi} of bivariate 
polynomials; the universal Grobner basis of the ideal / := ideal(P) is 

1 7 

U(I) = FU {iei — 3x? + 2xi, x% - hx\ + Ax 2 , x\ + -x\ — -x 2 , x 2 - x\\ 

bo 

and contains univariate polynomials in each of x\ and x%. Finding the roots of these two poly- 
nomials, we conclude that the set of zeros of P satisfies var(P) = var(J) C {0,1,2} x {0,1,4}. 
Substituting back to P, we find that var(P) = {(0,0), (1, 1), (2,4)}, thereby solving the system. 

The length of an ideal / in ¥[x] is the F-dimension of the quotient F[x]// and is finite if and 
only if the set of common zeros of / over the algebraic closure of F is finite. Let Hilb^ be the set 
of n-long ideals in ¥[x] = ¥[x\,... , xj\\ it can be embedded as an algebraic variety in a higher 
dimensional space and is referred to as the Hilbert scheme of n-long (f-variate ideals. One of 
the goals of this article is to provide a polynomial time algorithm for computing the universal 
Grobner basis of any ideal on the Hilbert scheme (see Section 4 for the complete formulation): 
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Theorem [4.2| Fix d. Then there is a polynomial time algorithm that computes the universal 
Grobner basis ti(I) of any ideal I £ Hilb^ using 0( n 2d+1 (log n)^ 2d ~ 1 ^ d ^ 1 ^ ) arithmetic operations. 

The computational complexity is measured in terms of the number of arithmetic operations over 
the underlying field F. Over the field of rational numbers, the algorithm is (strongly) polynomial 
time in the Turing computation model, but we do not dwell on the details here. 

One ingredient of our algorithm is an effective unified construction, for any ideal I £ Hilb^, of 
its state polyhedron S (I) whose vertices bijectively index the reduced Grobner bases of I. This is 
done in Section 2, where we introduce the basis polytope B(I) of any ideal I € Hilb^ and establish 
the following description of its state polyhedron (see Section 2 for the complete statement): 

Theorem ||]| The state polyhedron of / € Hilb^ is provided by S(I) := B(I) + Ml. 



As a corollary, we obtain the following polynomial upper bounds on the number of reduced 
Grobner bases and the size of the universal Grobner basis of any ideal on the Hilbert scheme: 

Corollary 2.5 For every fixed d, the following hold for any n-long ci-variate ideal / € Hilb^: 



,d-l 



(1) the number of distinct reduced Grobner bases of I is 0{n d + 1 

(2) the number of elements in the universal Grobner basis U(I) is 0(n^" " 1 ' d ( d + 1 )) 



2d-3+- 3d_1 



The cardinality of the set defining the basis polytope B(I) of / € Hilb d is typically exponential 



in n even for fixed d = 2 and so Theorem 2.4 does not lead directly to an efficient algorithm for 
constructing the state polyhedron. We overcome this difficulty by introducing, in Section 3, the 
Hilbert zonotope TC d . Proving that B(I) is a projection of a suitable matroid polytope, we show 
that Tt d is universal for the Hilbert Scheme in the following sense: 

Theorem |3.5| The Hilbert zonotope Ti d refines the state polyhedron S(I) of every / € Hilb^. 



Using Theorems and |3.5| we are able, in Section 4, to provide the aforementioned polyno- 
mial time algorithm for constructing the state polyhedron and universal Grobner basis IA{I) of 
any / E Hilb^. In particular, our results apply for the vanishing ideal of any point configuration, 
extending earlier results of for the generic case, and for lattice ideals studied earlier in []TTf . 

In Section 5 we interpret some of the notions and demonstrate some of the results discussed 
herein for the special classes of vanishing ideals of point configurations and of lattice ideals, the 
latter having some consequences for the so-called "group relaxation" of integer programming. 

We conclude with a brief discussion, in Section 6, of the embedding of the Hilbert scheme Hilb^ 
into the Grassmanian of n-dimensional subspaces of a vector space of dimension 0(n (log n) d ~ l ). 

2 The basis polytope and the state polyhedron 



A staircase is a set A C N rf of nonnegative integer vectors such that u < v <E \ (coordinatewise) 
implies u S A. Let ( N ) . . denote the finite set of n-element staircases in N d . For d = 2 the n- 

^ V n /stair 
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staircases are the Young diagrams of n. For a staircase A, let A := N d \A be its complement in N d 
and let min(A) be the unique finite set of coordinatewise minimal vectors in A. The n-staircases 
in N rf are in bijection with the monomial ideals in Hilb^ via I\ := ideallx 11 : v G min(A)}. 

Now fix any ideal / G Hilb^. An n-subset A C N rf is basic for I if the congruence classes 
modulo / of the monomials x v with v G A form a vector space basis for the quotient space F[x]/I, 
or equivalently, if the F- vector space linjx 1 ' : v G A} satisfies linjx 1 ' : v £ \}f^\I = {0}. If A is 
basic then the class [/] = / + / of any / G F[x] contains a unique representative in linlx 11 : v G A}; 
let [f]\ denote this unique polynomial satisfying [f]\ G linlx" : v G A} and / — [f]\ G L. 

A staircase A G ( \tair ls ™^ a ^ f or I ^ its monomial ideal I\ is the initial ideal in^(I) := 
ideal{in_ < (/) : / G /} of 7 under some monomial order -<;. If A is initial then it is also basic and 
the unique reduced Grobner basis of I under -< is the set G\(I) := {x u — [x u ]\ : u G min(A)} 
consisting of precisely |min(A)| polynomials. Let A(I) denote the set of initial staircases of /. 
We shall need the following two propositions on basic sets and initial staircases of an ideal. 

Proposition 2.1 Let -< be any monomial order and let A G A(J) be the initial staircase of I 
satisfying I\ = in^(I). Then for any vector u G N d \ A we have [x u ]\ G lin{x v : v G A, v -< u}. 

Proof. Let / := x u — [x u ]\ G I and let x v = in^(f). Then x v G in^(I) = I\ hence v G" A. Since 
x u is the only monomial in / with exponent not in A, it must be that v = u hence in^(f) = x u , 
so all monomials involved in [x u }\ are smaller than x u under -< as claimed. □ 



Proposition 2.2 Let -< be any monomial order and let A G A(J) be the initial staircase of 
L G HUb^ satisfying I\ = in^(I). Let fj, be any basic set of L. Write A = {Ai,...,A n } and 
\x = {{J,i, . . . , /i n } with Ai ~< ■ ■ ■ -< A n and /t*i [i n . Then A& ^ /ife /or each k = 1, . . . , n. 

Proof. Suppose indirectly /ife -< A^ for some k, and let U := lin{[x Ml ]A, • • • , [£ Mfc ]A/- Then 
A*i d: l^k -4 for alH < k hence [x^jx G lin{x Aj : j < k} (clearly if m G A and by Proposition 



2.1 if ^ G" A). So U C lin{x A J : j < /e} hence dim(f7) < k. But /x is basic so {[x^ 1 ]^, • • • , [^"Ia} 



is linearly independent hence so is {[x^Ja, • • • , [x^^Ja}- Thus, dim(C7) = k, a contradiction. □ 

Each vector w in the nonnegative orthant WL partially orders monomials x v by the value 
w ■ v. For any generic w this is a total order and hence a monomial order, and every monomial 
order arises that way from some generic w G The initial ideal of I under w is in w (L) := 
ideal{m U) (/) : / G 1} and is a monomial ideal if w is generic. Declare two non-negative vectors 
w and w' equivalent if in w (L) = in w >{L). The equivalence classes are relatively open convex 
cones forming a subdivision of called the Grobner fan of L (cf. H). A vector w lies in a full 
dimensional cone of the Grobner fan if and only if in w (L) is a monomial ideal. 

The (minimizing) normal cone of a face at a polyhedron P in M d is the relatively open cone 
of those vectors w G M. d uniquely minimized over P at that face. The collection of normal cones 
of all faces of P is called the normal fan of P. Let B be any polytope in M. d and let P := B + 
Then the normal fan of the polyhedron P forms a subdivision of and the vertices of P are 
precisely those vertices of B whose normal cone at B contains a strictly positive vector w G R+. 
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A polyhedron P in R d is the state polyhedron of the ideal I (cf. Q) if the Grobner fan of I 
equals the normal fan of P. This holds if and only if the set A(I) of initial staircases is in bijection 
with the vertex set of P, with w G uniquely minimized over P at its vertex corresponding to 
A G A(7) if and only if the initial monomial ideal of / under w satisfies in w (I) = I\. 

We now describe a construction of the state polyhedron for which this bijection is very nat- 
ural, in that the vertex corresponding to A G A(/) is simply its vector sum £A G N d . Let 
yd ■- \J( N n ) sta . r denote the union of all n-staircases in N rf . Given an ideal / G Hilb^, let 
r(I) := {A C V~ : A basic for /} denote the finite set of all n-subsets of V% basic for /. Since 

every initial staircase of / is basic and is contained in V d = U d \tair we nave A(/) — The 
following polytope will later enable the efficient computation of the state polyhedron. 

Definition 2.3 The basis polytope of I G Hilb^ is the convex hull of sums of basic sets of I in V d , 

B{I) := conv{^A : A GT(/)} C R d . 

The state polyhedron of any n-long d-variate ideal is provided by the following theorem. 
Theorem 2.4 The state polyhedron of I G Hilb^ is provided by S(I) := B(I)+M.i_. Furthermore: 

(1) its vertex set is {Y A : A G A(I)} and is in bijection with A(J) via the map A — ► Y A ; 

(2) a generic w G is minimized over S(I) and B(I) at ^ A with A G A(J) and I\ = in w (I) . 

Proof Let w G be any generic vector and let A G A (I) be the initial staircase of / 
satisfying in w (I) = I\. Consider any fj, G T(I). Writing A = {Ai, . . . , A n } and fi = . . . , fi n } 



with w ■ Ai < ■ • • < w ■ X n and w ■ m < ■ ■ ■ < w ■ u. n , we find by Proposition 2.2 that for each k 



we have w ■ A& <w ■ with equality if and only if A& = fik since w yields a total order. Thus, 

n n 

w ■ ^ = w ' ^ k — w ' ^ k = w ' ^ 

k=l k=l 

with equality if and only if A& = fik for all k. We find that for all n G T{I) other than A we have 
w-^2X < \x hence w is uniquely minimized over B{I) + at ^ A. This shows in particular 
that 5^ A is indeed a vertex of S(I) and B(I) and that Y2 f 1 Yl ^ f° r an y A 4 ^ other than 
A. Since any A G A(J) satisfies I\ = in w (I) for some generic w G K+, the map A — > Yl ^ 1S indeed 
a bijection from A(I) onto {Y A : A G A(/)}. Since generic vectors are dense in and a fan 
is determined by its maximal dimensional cones we find that the Grobner fan of I equals the 
normal fan of S(I) = B(I) + iH , showing that S(I) is indeed the state polyhedron of /. □ 

Theorem remains valid if we replace B(I) by the convex hull B = conv{ ^ A : A G T } 
with r ^ A(7) any collection of basic sets of I which contains all initial staircases of /. In 
particular, it holds with F = A(I) and with T = {A G (^) . . : A basic for /}. However, these 
choices do not lend themselves to efficient algorithmic construction of the state polyhedron and 
the universal Grobner basis. As we shall see, we need to take T to be the collection of all basic 
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sets of I contained in some subset V C N rf . The smallest such set yielding T which contains all 
n-staircases in N rf is their union V^, leading to our choice of T(I) and £>(/). 

Let IA{I) denote the universal Grobner basis of I defined as the union of all reduced Grobner 
bases of I and hence simultaneously providing a Grobner basis for / under any monomial order. 



Theorem 2.4 yields the following polynomial upper bounds on the number of distinct reduced 
Grobner bases and the size of the universal Grobner basis of any ideal in Hilb^. The proof is 



similar to that of fllO[ ] for the case of vanishing ideals of point configurations. 

Corollary 2.5 For every fixed d, the following hold for any n-long d-variate ideal I G Hilb^: 

(1) the number of distinct reduced Grobner bases of I is 0(n d + 1 ) ; 

(2) the number of elements in the universal Grobner basis U\I) is 0{n <*(<*+i) ) . 

Proof. Consider any A £ ( H n ) . . Then, for % = 1, . . . , d, each v G A satisfies vi < n hence the 
ith coordinate value of ^ A is less than n 2 . Thus, the lattice polytope P = conv{^ A : A € A(/)} 

„ , d — 1 d—1 . . 

is contained in the cube [0, n ] and hence has 0(vol(P) d+1 ) = 0(n d + 1 ) vertices (cf. By 



Theorem 2.4, the state polyhedron S(I) = P + and P have the same vertex set, giving (1). 

Next, for any A € ) . . , the size of the set min(A) of minimal elements not in A is 0{n~3~) 
(cf. fH). Since the reduced Grobner basis G\ of I corresponding to A has | min(A)| elements, the 
product of the bound on any | min(A)| and the bound just established in (1) on the number of 
reduced Grobner basis yields the bound in (2) on the size of their union U{I). □ 



3 The Hilbert zonotope and its universality 



While the number of initial staircases of J G Hilb^ is polynomial in n for any fixed d by Corollary 



2.5, the cardinality of T(I) is typically exponential in n even for d = 2; for instance, if I is the 
bivariate vanishing ideal of any generic n points in the plane over an infinite field (see discussion 
of point configurations in Section 5), then all n-subsets of V£ are basic, and hence so are all 
n-staircases in N 2 which are in bijection with number partitions of n. Thus, it is not possible to 
filter A(J) out of T(I) and construct B(I) or S(I) directly in polynomial time. 

To overcome this we now introduce, for each pair of positive integers d and n, the Hilbert 
zonotope TC^. As we shall see, it is universal for the Hilbert scheme Hilb^ in that it provides a 
refinement of the basis polytope B(I) and the state polyhedron S(I) of every ideal I G Hilb^. 

Recall that V? = I J d h ■ is the union of all n-staircases in N d . Call an element of the 

n w V n /stair 

symmetrization — V£ = {u — v : u,v G V£ } primitive if it is not a nonnegative integer 
multiple of another element of V!~ — V£- For n > 2 let be the set of primitive elements of 
Vn — V£, and for n = 1 let Df := ±{ei, . . . , e^}. We make the following fundamental definition. 
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Definition 3.1 The Hilbert zonotope TC d is the following Minkowski sum of line segments, 

U d n := £>l].i; C R d . 

veDi 

Note that = —D d is centrally symmetric, hence so is 7i d which equivalently can be defined 
as 1, 1] • v by summing over only one of each pair {— v,v} of antipodal primitive elements. 

The Hilbert zonotope is well behaved in the following sense. 

Proposition 3.2 Fix any d. Then the number of vertices of the Hilbert zonotope TC d satisfies 
Q^ n 2(d—i) ^ Q g n ^2(d—l) ^. f ur t} ier) i n polynomial time using that many arithmetic operations, all 
its vertices can be listed, each h along with a vector w(h) uniquely minimized over TL d at h. 

Proof. First note that V d = (J (^) stair is given by V d = {v G N d : IliLiC^ + x ) < n Y 
indeed, if v G N d lies in some n-staircase A then the entire box {u G N d : u < v} is con- 
tained in A hence riiLi^i + 1) = ll n G N d : n < w}| < |A| = n; conversely, if v G N rf sat- 
isfies nf=i( u i + -0 — n then the box above can be augmented with multiples of unit vectors 
to an n-staircase A containing v. Now the cardinality of {v G N rf : nf=i( w « + 1) — n l obeys 
the bound 0(n (log n) rf_1 ) (cf. Jl^]). Thus, the Hilbert zonotope TC d is the Minkowski sum 
of N := %\D*\ = 0(\V d \ 2 ) = 0(n 2 (log n) 2 ^ 1 )) line segments and therefore (cf. Q, g) has 
0(N d ~ l ) vertices which can all be enumerated, each h along with a vector w(h) uniquely mini- 
mized at h, using 0(N d ~ 1 ) arithmetic operations, giving the claimed bounds. □ 



Example 3.3 We compute the Hilbert zonotope Ti^ for d = 2, n = 3. Using Proposition 3^ we 
find the set H of vertices of TL 2 which shows that it is a (centrally symmetric) 10-gon, as well as 
the corresponding set W of the linear functionals uniquely minimized at its various vertices: 
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We need to recall a few facts about matroids and matroid polytopes. Let M = (V, T) be a 
matroid over a finite set V with collection of bases T C 2 V . Its matroid polytope is defined as the 
convex hull B{M) := convjl^ : BeT}c R , where 1b denotes the incidence vector of B C V, 
that is, the {0, 1} vector in M v with support B. This is a well known object of importance in 
combinatorial optimization. Below, let e v G R v denote the unit vector indexed by v G V . 



Proposition 3.4 Every 1-face of the matroid polytope is equal to e u — e v for some u,v 6 V. 
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Proof. Consider any pair A, B G T of bases such that [1a,1_b] is an edge (that is, a 1-face) 
of B(M), and let w G R^ be a linear functional uniquely maximized over B(M) at that edge. If 
A \ B = {u} is a singleton then B \ A = {v} is a singleton as well in which case 1a — 1b = e u — e v 
and we are done. Suppose then, indirectly, that it is not, and pick an element u in the sym- 
metric difference AAB := (A \ B) U {B \ A) of A and B of minimum value w u . Without loss 
of generality assume u G A \ B. Then there is a v G B \ A such that C := A \ {u} U {v} 
is a basis of M. Since |ylA£>| > 2, C is neither ^4 nor 5. By the choice of u, this basis satis- 
fies w ■ lc =io- 1a — "U^m + > u> • 1a, hence also a maximizer of w over B(M), a contradiction. □ 

A polyhedron P is a refinement of a polyhedron Q if the normal fan of P is a refinement of 
that of Q, that is, the closure of each normal cone of Q is the union of closures of normal cones 
of P. The significance of the Hilbert zonotope is now demonstrated by the following theorem. 

Theorem 3.5 The Hilbert zonotope TL d is a refinement of both the basis polytope B(I) and the 
state polyhedron S(I) of every member I of the Hilbert scheme Hilb d of n -long d-variate ideals. 

Proof. Consider any / G Hilb^. Let M := (V^,T(I)) be the matroid over with collection 
of bases r(/), which is the restriction to V£ of the infinite matroid over N d of all basic sets of I. 
Let B(M) := conv{l A : A G r<7)} C R v * be the matroid polytope of M and let 

vr : R v " — ► R d : e v i-> v 

be the natural projection sending the unit vector e v corresponding to v G to the vector 
v eN d cR d . Then for each A G T(J) we have tt(1\) = Y) A hence the basis polytope of / 

B{I) = conv{^A : A G T(I) } = conv{vr(l A ) : A G T(J) } = tt(B(M)) 

is a projection of the matroid polytope. Thus, each edge of B(I) is the projection of some edge 
of B(M) hence, by Proposition |3.4| , is equal to 7r(e n — e v ) = u — v for some pair u,v G V£ and 
therefore parallel to some element in D d . Thus, the Hilbert zonotope H d = ^2 v£ £>d[0, 1] • v is 
the Minkowski sum of a set of segments containing all edge directions of B(I) and therefore its 
normal fan is a refinement of the normal fan of B(I). 

Next, consider any face F of S(I). Then F is also a face of B(I) C S(I) and hence the closure 
of its normal cone at B(I) is a union of closures of normal cones of Tt d . But D d contains all unit 
vectors which implies that each normal cone of Ti. d is contained in the interior of some orthant. 
Thus, the closure of the normal cone of F at S{I) is the union of the closures of those normal 
cones of 7i d contained in the normal cone of F at B(I) which lie in the nonnegative orthant R^. □ 

Let us call the coarsest common refinement of the state polyhedra of all ideals / on the 
Hilbert scheme the Hilbert polytope. As pointed out to us by D. Bayer, the Hilbert zonotope may 
in general be finer than the Hilbert polytope; in particular, it may have more vertices. However, 
the Hilbert zonotope allows efficient algorithmic treatment while the Hilbert polytope might not. 
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4 Computing the state polyhedron and universal Grobner basis 

We now use the Hilbert zonotope to efficiently compute the set of initial staircases, the state 
polyhedron and the universal Grobner basis of any ideal on the Hilbert scheme. 

Let U% := {u + et : «£ V~, < i < d} with eo := and ej the ith unit vector for 1 < i < d. 
Then U% contains, along with every ra-staircase A, all vectors u + with u any vector in A and 
any unit vector, hence also min(A). On the other hand, Uf° C V^ n and hence, for fixed d, obeys 
the same upper bound 0(n (log n) d_1 ) on its cardinality as does V^. 

We assume that / G Hilb^ is presented by its Grobner basis G\ := {x u — [x u ]\ : u G min(A)} 
under some monomial order -<. Such a presentation, say with respect to the degree reverse 
lexicographic order, is known to be efficiently computable from any generating set - see discussion 
at the end of this section. Given such an ideal we will need the representative [x u ]\ for every u 
in the set U d . We include the short proof of the following adaptation of || Propoistion 3.1]. 

Proposition 4.1 Fix any d. Then, given any ideal I = ideal(G\) G Hilb^, the representatives 
[x u ]\ = ^2 v( z\ a v ^ u -x v of all u G U% can be computed using 0(n 3 (log n)^ 1 ) arithmetic operations. 

Proof. Compute [x u ]\ for the elements u S Uf t in the order -< as follows. If u € A then 
[x u ]\ = x u . If u E min(A) then x u — [x u ]\ is in the given reduced Grobner basis G\ from 
which [x u ]\ can be recovered. Otherwise s := u — G U% \ A for some i. By Proposition 2.1 
[x s ] x G lin{x* : t G T} with T := {t G A : t -< s}. Now s -< u and t + ej -< s + e« = u for all t G T, 
hence [x s ]\ = J2teT a t,s xt an d [^* +<3i ]a = J2ve\ a v,t+e t x v for all t G T are already available. Thus 

[x u }\ = [x s+e % = ^2a t , s [x t+e >} x = ^2a ttS ^2a V)t+ei x v = ^(^at^t+e^ 

teT tgT vex vex teT 

is now obtained using 0(n 2 ) arithmetic operations. Since \U^\ = 0(n (log n) d_1 ) we are done. □ 

We are finally in position to provide the efficient procedure for computing the set of initial 
staircases, the state polyhedron and the universal Grobner basis of any ideal on the Hilbert 
scheme. Theorem |3.5| enables us to bypass the difficulty caused by the exponential size of r(I). 



Theorem 4.2 For every fixed d there is a polynomial time algorithm that, given any n-long d- 
variate ideal I G Hilb^, computes its set A(I) of initial staircases, its state polyhedron S(I), and 
its universal Grobner basis using 0(n 2d+1 (log n)^ -1 ^ -1 )) arithmetic operations. 

Proof. First enumerate, as in Proposition and hence within the claimed complexity bound, 
the 0(n 2 ( rf_1 ) (log n) 2 ^ 1 ) ) vertices of Ti.^, each vertex h along with a vector w(h) uniquely 
minimized over TC d at h. We claim that any (coordinatewise) positive w(h) on the list is uniquely 
minimized over B{I) at ^ (i for some initial staircase u G A(/) and, conversely, for every fi G A(J) 
there is some (possibly many) positive w(h) on the list uniquely minimized over B(I) at ^ \x. 
Consider first any positive w(h) on the list. Since Ti^ refines B(I) by Theorem |3.5| , the vector 
w{h) and hence a positive generic perturbation of it lie in the normal cone of some vertex of £>(/). 
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By Theorem 2.4, this vertex is for some initial staircase fi G A(I). Conversely, consider any 
initial staircase fi G A(J). Then, by Theorem |2.4| , ^ /i is a vertex of 5(7) hence its normal cone 
is contained in the interior of the nonnegative orthant. Since TC^ refines S(I) by Theorem 3.5, 
this normal cone contains the normal cone of some (possibly many) vertex h of TC^- The vector 
w(h) of that vertex h on the list must then be positive. 

We proceed to show that for each positive w(h) on the list we can efficiently compute the 
minimizing vertex ^ /i of B{I), the initial staircase //, and the reduced Grobner basis G^. 

First, compute, as in Proposition 4.1, using the given presentation / = \deal{G\) of / by some 
reduced Grobner basis, the representative [x u \\ = X^eA a v,u • x v of every vector u G U%. 

Now pick any vertex h of TL^ with w := w(h) positive. Use the following greedy algorithm to 
find a basic set /j, = . . . , G r(J) : for i = 1, . . . , n pick to be an element in "l^f C U% of 
minimal value w ■ \i{ with the property that {[x^Ja ■ j <i} is linearly independent. Then for any 
other basic set v G writing v = {v\, . . . , u n } with w-u\ < ■ ■ ■ < w- v n , we have w • ^ < w ■ 
for all k hence w ■ ^ jj, < w ■ ^ v. This shows that ^ /i minimizes w; over B(I) hence, as shown in 
the first paragraph of this proof, is the unique minimizing vertex of the positive w(h) = w over 
13(1), and ji G A(I) is the corresponding initial staircase. 

The determination of \x by the greedy algorithm can be efficiently implemented by performing 
Gaussian elimination on the fly as follows. Totally order compatibly with w so that if u 
precedes v then w ■ u < w ■ v and let A = {Ai, . . . , A n } be an arbitrary labeling of A. Thus, the 
coefficients a„ jU for all u G and v G A now form an n x matrix A over F. For each u G 
let A u G F n denote the column of that matrix corresponding to u, with Af := a\ itU . Now, in 
the greedy algorithm, for i = 1, . . . , n, pick fii to be the first element in C U% whose column 
A^ contains a nonzero coordinate A^ for some k > i. Apply suitable row operations to A so as 
to transform A^ to the unit vector while maintaining = ej for all j < i. This consumes 
0(n ■ \Un\) arithmetic operations per /ij, totaling to 0(n 2 ■ \U^\) operations. 

Moreover, the updated matrix A at the end of this iterated process is the matrix of coordi- 
nates in the new basis ^: for each u G U^, the representative [x n ]^ is simply read off from the 
corresponding column as [x u ]^ = ^7=1 Afx^. Now, the set min(/2) consists precisely of those 
u G U% \ A 4 with the property that for each i = 1, . . . , d either ui = or u — G fx, hence can be 
quickly filtered out of U%- So, within the same complexity bound of 0(n 2 ■ \U%\) operations we 
obtain the reduced Grobner basis G M := {x u — [x u ]^ : u G min(jl)} corresponding to [i. 

Summarizing, for each positive w(h) on the list of vertices h of TC^, we can find the ini- 
tial staircase fi G A(/) for which ^ /j, is the unique minimizer of w(h) over B(I) and the 
corresponding reduced Grobner basis using 0(n 2 ■ \U^\) operations. Therefore, the entire 
set A(I) of initial staircases, the state polyhedron S(I) = conv{^ [i : fi G A(/)} + R^, and 
the universal Grobner basis U (I) = U At gA(/) °f the given ideal / can be produced using 
Q^ n 2(d-i) ^ Q g n ^2(d-i) 2 . n 2 . n n )d-i) ) arithmetic operations as claimed. □ 



Example 4.3 We compute the universal Grobner basis of the ideal / = ideal(GA) G Hilb^ for 
d = 2, n = 3 over F = R, presented by G\ = {x\ — 3xf + 3xi — 1, x% — x\ + 1} which is 
its reduced Grobner basis under the lexicographic order with X2 > x\, with A the staircase 
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A = {00, 10, 20}. First, we obtain the set of w £ via the Hilbert zonotope as in Example 
|3.3| , which is W+ := {31, 32, 23, 13}. Next, we compute the 3 x 10 matrix A of coefficients of the 
representatives [x u ]\ = a\ )U ■ 1 + a^u • x\ + a^,^ ■ x\ for all u G C/f , whose rows and columns are 
indexed by A and C/f respectively, and obtain 



1 1 

-2 -2 
1 1 

compatibly with w and suitably 
permute the columns of A. Apply the greedy algorithm and find the new initial staircase u = 
{00,01,02}. Next apply suitable row operations to A to make the transformation to the new 
basis and obtain the following updated matrix, with rows and columns suitably re-labeled, 





00 


10 


20 


30 


01 11 21 


00 


( 1 








1 - 


-1 1 


A = 10 





1 





-3 


1 -1 -3 - 


20 







1 


3 


1 2 


some vector in 


w+, 


say 


w : 


= 32 


Reorder £/| 



00 


01 


10 


02 


11 


03 


20 


12 


21 


30 


1 





1 




















1 





1 


1 





1 





2 





1 


3 











1 


1 





1 


1 


2 


3 



A = 

Now, we have min(/i) = {10, 03} and so the new reduced Grobner basis is read off from the third 
and sixth columns of the new matrix to be = {x\ — X2 — 1, x|}. 

Repeating this for each of the other three vectors in W+, we keep getting either A or jj,. 
We conclude that here A(7) = {A, fi} consists of two staircases only, the state polyhedron is 
S(I) = conv{^ A, ^2 ^} + ^+ = [30, 03] + and the universal Grobner basis is 

U(l) = GxUG^ = {x?-3xf + 3xi-l, x 2 -x 1 + l, x x -x 2 -l, x 3 2 } . 

It is known (cf. || and references therein) that the reduced Grobner basis under the degree 
reverse lexicographic order of any ideal / = ideal(-F) G Hilb^ presented by any set of generators 
can be computed in time 0{D d2 ) where D is the maximal degree of any generator / G F. 



Thus, Theorem 4.2 also implies, for any fixed d, an efficient algorithm for computing the state 



polyhedron and the universal Grobner basis of any ideal of the Hilbert scheme Hilb^ presented 
by any set of generators using a number of arithmetic operations polynomial in n and D. 



We conclude by pointing out that Theorems [3J} and |4.2| yield the practical outcome that, 
for each d and n, a list W% of positive w(h) of vertices h of TC d can be computed once and for 
all, providing a universal set of monomial orders for the Hilbert scheme Hilb^ that allows the 
efficient computation of the universal Grobner basis hi (I) of any given n-long d-variate ideal /. 

5 Examples: point configurations and lattice ideals 

We now interpret some of the notions and results discussed above for several special classes of 
ideals of the Hilbert scheme Hilb^. We start with the simple class of monomial ideals. 
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Monomial Ideals. Recall that the monomial ideals in Hilb^ are in bijection with the n-staircases 
in N d via I\ = idealjx" : u G min(A)}. Consider any such ideal I\. Then T(I\) = {A} is a 
singleton: indeed, if u is any other n-subset of N d and u G u \ A then x u G linjx*' : !) £ fl}f)Ij 
hence \i is not basic. For any monomial order -< we have in^(I\) = I\ hence A(J\) = {A} as well. 
Thus, B(I\) = {^2 A} is a single point and the state polyhedron is S(I\) = {^2 A} + For 
every u G N d \ A we have G 7> hence [x u ]\ = 0. Thus, the universal Grobner basis of 7 equals 
the unique reduced Grobner basis G\ and are both given by U(I\) = G\ = {x u : u G min(A)}. 

Point Configurations. The vanishing ideal Iq '■= {/ G ¥[x] : f(c\) = ■ ■ ■ = f{c n ) = 0} of a 
configuration C = {c\, . . . , c n } of n distinct points in affine space ¥ d is a radical ideal of length 
dim(F[x]/7c<) < n. Assume throughout this example that F is infinite, which implies that the 
length of Ic is exactly n and hence Ic G Hilb^. For A = {Ai, . . . , X n } C N d let 



/c Al 
X 



\ \ \ 



where c^ 3 = c^ 3 ' k 



k=l 



\ r^ 1 r^' 1 ■ ■ ■ r^ n I 

Then A is basic for Ic if and only if C A is nonsingular: indeed, the vector (/(ci), . . . , f(c n )) G F n 
of evaluations of a polynomial / = Y^2=i a i xXi ^ lin{x Al , . . . ,x Xn } at the points of C is provided 
by C A ■ A and is the zero vector in F n if and only if / lies in 7; thus, lin{x Al , . . . , x An } n Ic = {0} 
if and only if C A is nonsingular. Thus, T(Ic) = {X C VJ~ : det(C A ) / 0}; the basis polytope is 
B(I C ) = conv(]P A : A G r(7 c )}; the state polyhedron is S(I C ) = B{I C ) + M+; and the set of 
initial staircases is A (7c) = {A G r(7c) : Yl A vertex of S (Ic) }■ For every sufficiently generic 
configuration, say one satisfying det(C A ) ^ for all A C V d , the state polyhedron S(Ic) and 
the set of initial staircases A(7c) coincide, respectively, with the corner cut polyhedron P d and 
the set ) of n-element corner cuts in N d introduced and studied in |l(| . If -< is a monomial 
order and A G A(7c?) is the initial staircase with in^(Ic) = I\ then for every u G N rf we have 
[x u ] x = x u - det( ({x} U C)^"> UA ) ) • det -1 (C A ). So the reduced Grobner basis of I c under -< is 

G A = {det( ({x} U C) ({n}uA) ) • det _1 (C7 A ) : u G min(A) } 

and the universal Grobner basis U(Jc) = U\eA(i c )G\ is efficiently computable by our methods. 
This extends the results of [pi]] where an efficient construction oilA(Ic) for generic configurations 
only was provided, based on the separable-partitions methods of |2], [?J • 

Lattice Ideals. The binomial ideal of an integer lattice L C Z rf is 7l := idealjx 1 ^ — x v : v G 7} 
where v + ,v~ G N d denote the nonnegative and nonpositive parts of v G N rf with i> = t> + — v~ . 
Assume L is full dimensional with determinant det(L) = n implying that II has length n hence 
II G Hilb^. An n-subset A = {Ai, . . . , A n } is basic for It, if and only if it is a set of distinct 
representatives of the congruence classes of Z d modulo L, that is, if and only if for all i ^ j 
we have Aj — Xj 7: this follows from the fact that u — v G L if and only if x u — x v G 7l 
(cf. 0). Thus, T(7 L ) = {A C V* : i ^ j implies Ai - Xj L}; the basis polytope is B(I L ) = 
conv(^A : A G r(7^)}; the state polyhedron is <5(7l) = B(It) + and the set of initial 
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staircases is A(Jt) = {A € T(Il) : ^ A vertex of S(Il) }. If -< is a monomial order and A € A(ii) 
is the initial staircase with in^(Iy) = -T\ then for every u G N d we have [x"]^ = x" A with u\ 
the unique representative in A with u — u\ G L. So the reduced Grobner basis of under -< is 
G\ = {x u — x Ux : u G min(A)} and the universal Grobner basis U{Il) = Uxe\(i L )G\ is efficiently 
computable by our methods. Here the universal Grobner basis consists of binomials only, and 
the set of integer vectors T := {u — v : x u — x v € U{Il)} which can be read off at once from 
U{l£) is a universal test set for the lattice minimization problem that, given any x G N d and any 
w G asks for an x* G N rf satisfying x — x* G L and minimizing the value w ■ x. 



6 On the coordinatization of the Hilbert scheme 

We conclude with a brief discussion of the embedding of the Hilbert scheme Hilb^ into the 
Grassmanian of n-dimensional subspaces of a vector space of dimension 0{n (log n) d_1 ). 

Recall V d = u(^ ) stair is the union of n-staircases and U d = {v + &i : v G V d and < i < d} 
with eo = and the i-th unit vector in M. d for 1 < % < d. Throughout this section, assume that 
U d is totally ordered. Let S := ¥[x u ...,x d ], let Syd := lin{x u : u G U%} C S be the F-linear 
span of the monomials x u with u G U d , and for each ideal / G Hilb^ let lyd := I D Syd. 

Proposition 6.1 The dimension of the vector subspace lyd of any ideal I G Hilb^ isp := \U d \—n. 

Proof. Consider the map (ft : Syd — > (S/I) such that / = Ylueu d a uX u >—>/ + /. Then 
ker(^) = {/ G Syd : / G 1} = /[/d. Hence the map ip : (Syd/Iyd) — > (S/I) such that 
/ + lyd h-> / + / is injective. Consider any element 5 + / of 5/7 and let [^]a G Sj/d be the normal 
form of g with respect to the reduced Grobner basis of / corresponding to some initial staircase 
A of I. Then ip([g]\ + lyd) = [g]\ + 1 = g + 1 hence if) is surjective. Therefore, S/I and Syd/Iyd 
are isomorphic F-spaces hence dim(Syd/Iyd) = dim(5/I) = n and dim(Iyd) = \U d \ — n = p. □ 

Proposition 6.2 Any ideal I G HUb^ is uniquely determined by its vector subspace lyd. 

Proof. By definition of U d , every reduced Grobner basis of / lies in lyd, implying that the ele- 
ments of lyd generate / as an ideal. Thus, for any /, J G Hilb^, / = J if and only if lyd = Jyd. □ 



Proposition |0] shows that lyd is a p-dimensional subspace of Syd, i.e., a point on the Grass- 
manian Gv(Syd,p) of p-dimensional subspaces of the |[/^ (-dimensional vector space Syd. This 



implies that lyd, and therefore, by Proposition p.2| , also /, inherits standard Plucker coordinates 



from the Grassmanian, leading to an embedding of the Hilbert scheme in projective space, as 
follows. Let {fi) p i=l be any ordered vector space basis for lyd with fi = ^2 uG yd rrii^ u x u , and for a 
p-subset v C N d let M u := (m i; 

u)i<i<p,ueu be the corresponding p x p submatrix of coefficients. 
The Hilbert scheme is then embedded into projective space by 



mibi — ► PC?) : 1^(1) := { det(M„) : v G j 
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it is well known that the Pliicker point (I) is independent of the choice of basis (/j) of Ijjd. 

We proceed with the dual embedding into the Grassmanian Gr (Sjjd,n) of n-dimensional 
subspaces of Sjjd. The space Sjjd is endowed with the standard monomial basis {x u : u € U^} 
and corresponding standard inner product {x u ,x v } = 5 U:V . Let I^ d C Sjjd be the orthogonal 
complement of Ijjd in S V d, which is isomorphic to Sud/Ijjd. Let (/ii)f =1 be any ordered vector 
space basis for I^ d with hi = J2 u eu d a i,uX u , and for an n-subset A C U% let A\ := (aj ;U )i<j< njUg A 
be the corresponding n x n submatrix of coefficients. The dual embedding is then given by 

(■) ± : Hilb£— ►p('$) : := jdet(A A ) :AG^" 

Next we explain how to actually compute the (dual) Pliicker coordinates; note, though, that 
the number of coordinates is Q^C ^™) ) hence exponential, so this computation can not be 
carried out in polynomial time even for d = 2. Let I = ideal(GA) £ Hilb^ be presented by 
its Grobner basis corresponding to some initial staircase A = {Ai, . . . A ra } under some monomial 



order -<. Compute as in Proposition LI the representative [x u ]\ = Y17=i a i,uX Xt of every u £ U^, 
and for i = 1, . . . ,n let hi := Yl u eu d a i,uX u - Then (hi)f =1 is an ordered basis of I^ d . With £/~ 
assumed to be totally ordered, we get the n x \U%\ matrix A = (aj )U ) of coefficients of this basis; 

the Pliicker coordinates are then read off from the minors of A as {I) ± = {det(j4^) : \i £ { U ")}- 
In particular, the Pliicker coordinate (I)j^ := det(j4^) is nonzero if and only if \x is basic for /. 

For the following classes of ideals, the Pliicker coordinates have a natural simple form. 

Monomial Ideals. Let I\ = \dea\.{x v : v £ min(A)} be the monomial ideal in Hilb^ correspond- 
ing to an n-staircase A. Then (see Section 5) A is the only basic set of I\. Therefore the only 

i I u d l ~\ 

nonzero Pliicker coordinate is and so (I) = e\ is the unit vector in projective space « > . 

Point Configurations. Let Iq be the vanishing ideal of a configuration C = {c±, . . . , c n } CF d 
with F infinite. For i = 1, . . . , n let hi := Ylu£U d c 7 xU with c" = Y\t=i c 7k- ^ or an y polynomial 
/ = J2ueu d m uX u £ Sjjd, its inner product with hi satisfies (/, hi) = YlueU d m utf = f{ci), hence 
is zero if and only if / vanishes on q. Thus, / is orthogonal to lin{/ii, . . . , h n } if and only if it 
vanishes on C, or equivalently, / £ Ic- This shows that (h{)f =1 is an ordered basis of (Ic)jjn. 

Let A = (cLi. u ) be the n x \U^\ matrix whose rows are indexed by the points c\, . . . ,c n in C 
and columns by u £ U%, with aj jU := c". Then A is the matrix of coefficients of the basis (/ij)™ =1 
and hence the Pliicker coordinates (Ic) of Ic can be read off from the minors of A. 
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